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fl, ■ Abstract: 

The longstanding problem of explaining the observed polarization of A hyperons 
inclusively produced in the high energy collisions of unpolarized hadrons is tackled by 
considering spin and k± dependent quark fragmentation functions. The data on A's 
and A's produced in p — N processes are used to determine simple phenomenological 
expressions for these new " 'polarizing fragmentation functions" , which describe the 
experiments remarkably well. 



1. Introduction 



It is well known since a long time that A hyperons produced with xf ^ 0.2 and 
Pt ^ 1 GeV/c in the collision of two unpolarized hadrons, AB — > A^X, are polarized 
perpendicularly to the production plane, as allowed by parity invariance; a huge 
amount of experimental information, for a wide energy range of the unpolarized 
beams, is available on such single spin asymmetries P]: 



da 



da 



AB—*A^X 



da AB ^ x + da AB ~ 



>AlX 



(1) 



Similar effects have been observed for several other hyperons, but we shall consider 
here only A's and A's. 

Despite the wealth of data and the many years they have been known, no con- 
vincing theoretical explanation or understanding of the phenomenon exist [[| [J. The 
perturbative QCD dynamics forbids any sizeable single spin asymmetry at the par- 
tonic level ; the polarization of hyperons resulting from the strong interaction of 
unpolarized hadrons must then originate from nonperturbative features, presumably 
in the hadronization process. A number of models attempting some understanding 
of the data in this perspective @-|| only achieve partial explanations. 

In the last years other single spin asymmetries observed in p^p — > ttX reactions 

]TTJ-|f20|; a phenomenological de- 
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have attracted a lot of theoretical activity 
scription of such asymmetries appears now possible with the introduction of new 
distribution |21|, [□], [L4], |22[ and/or fragmentation [JTJ], |TJ, |23|] functions which are 
spin and k± dependent; k± denotes either the transverse momentum of a quark 
inside a nucleon or of a hadron with respect to the fragmenting quark. 

that is, the azimuthal angle 



In particular the effect first discussed by Collins [12 



dependence of the number of hadrons produced in the fragmentation of a transversely 
polarized quark - has been recently observed [|4], 25] ; were such results confirmed 
the role of these new fragmentation functions would be of great phenomenological 
importance. 

We consider here an effect similar to that suggested by Collins, namely a spin 
and k^ dependence in the fragmentation of an unpolarized quark into a polarized 
hadron: a function describing this mechanism was first introduced in Ref. |23| and 
denoted by D^ T . This function is introduced in a frame defined by two light-like 
four- vectors n + and n_, satisfying n + ■ n_ = 1, and by the plane transverse to 
them. The four-momentum P of the outgoing hadron - a A hyperon in the present 
investigation - is in the n_ direction (up to a mass term correction). The function 
D^ T is then defined as (displayed in the n + ■ A = gauge) 



4 k Ti s Tj D ^ (zk 



x Tr (0\iP{y)\P, S T ; X)(P, S T] XF0(O) T -|O> 
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where the final state depends on the transverse part (St) of the spin vector S of the 
produced A only, i.e. one should interpret it as \P, St', X) = (\P, S = St', X) — \P, S = 
—St', X))/2, such that the contribution from unpolarized fragmentation cancels out. 
Furthermore, k± = \k±\ is the modulus of the transverse momentum of the hadron 
in a frame where the fragmenting quark has no transverse momentum. More details 
on this type of definition of fragmentation (or decay) functions can be found in Refs. 



MUM- 

In the notations of Ref. [IS] a similar function is defined as: 



A N D h y a (z, k ± ) = D h y a (z, k±) - D h i/ a (z, fcj.) (3) 
= D h y a (z,k±) -D h y a (z, -k±) , 

and denotes the difference between the density numbers D h y a (z, kj_) and D h i/ a (z, 
kj_) of spin 1/2 hadrons h, with longitudinal momentum fraction z, transverse mo- 
mentum k± and transverse polarization f or j, inside a jet originated by the frag- 
mentation of an unpolarized parton a. From the above definition it is clear that the 
k± integral of the function vanishes. 

The exact relation between D^ T and A N D h r i a is given by (notice that also D^ T 
should have labels h and a which are often omitted): 

A N D h y a (z, k ± ) = sin <P Di T {z, k ± ) , (4) 

ZM h 

where <fi is the angle between k± and the transverse polarization vector of the hadron, 
which shows that the function A N D h y a (z, k±) vanishes in case the transverse mo- 
mentum and transverse spin are parallel. 

In the sequel we shall refer to A N D h ^ i a and D\ T as "polarizing fragmentation 
function^' . 

In analogy to Collins' suggestion for the fragmentation of a transversely polarized 



quark we write |12L |27j| : 



D h y q (z, k ± ) = \ D h/q (z, k x ) + \ & N D h y q (z, k ± ) P \ ' {Pq * k±) (5) 

L L \p q X ACj_| 

for an unpolarized quark with momentum p q which fragments into a spin 1/2 hadron 
h with momentum p h = zp q + and polarization vector along the f = Ph direc- 
tion; Dh/q( z ,k 1 _) = D h y q (z, k±) + D h i/ q (z, kj_) is the kj_ dependent unpolarized 
fragmentation function. Notice that Ph ■ (p q x k±) = p q ■ (k±_ x Ph) ~ sin0. 

A QCD factorization theorem gives for the high energy and large pt process 
AB — > A^X, at leading twist with collinear parton configurations: 



da AB^x = j2 f a/A ( Xa ) g, f b/B ( Xb ) g, da ab ^ d ® D A y c (z) (6) 

a,b,c,d 
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and 

da AB^Mx = f a/A ( Xa ) <g, f b/B ( Xb ) <g, da" D Al/c {z) . (7) 

a,b,c,d 

Here and in the sequel we shall fix the scattering plane as the x-z plane, with hadron 
A moving along +z and the detected A produced in the first x-z quadrant; the |, 
I directions are then respectively +y and — y. 

In the absence of intrinsic k± (or rather when integrated over) the fragmentation 
functions _D a t i c {z) = J d 2 k±D A ^ / q (z, k±) (or D A i i c {z)) cannot depend on the hadron 
polarization, so that one has da^ = da\ which implies Pa = 0. 

However, by taking into account intrinsic k^ in the hadronization process, and 



assuming that the factorization theorem holds also when fc^'s are included |12 |, one 
has, using Eq. (|5|) instead of D A / C (z) in Eqs. ([I]), (^) and (|7|): 

E A d 3 a AB ^ AX „ ^ fdx a dx b dz 



d 3 P/ 



P\ = J T Z d * k± fa/A{Xa) fb/B(Xb) 



ab^cd 



x s6(s + i + u)— j — (x a , x b , k±) A^D^Jz, k±) (8) 
dt 



where s, t and u are the Mandelstam variables for the elementary process: for 
collinear configurations s = x a x b s, t = x a t/z and u = x b u/z and the modifica- 
tions due to intrinsic k± will be taken into account in the numerical evaluations. 
E A d 3 a AB ^ AX /d 3 p A is the unpolarized cross-section: 

E A d 3 a AB ~* AX ^ f dx a dx b dz 2 



d 3 ps 



Ef dX a U,X b UZ n I r ( \ r ( \ 

J ~ 2 d k ± f a /A{Xa) fb/B{Xb) 

u,,b,c,d 

^Q-ab-^cd 

x §S(s + i + u) - 2 — {x a ,x b ,k ± ) D A/c (z,k±) . (9) 

dt 



In Eq. (H) k± is considered only where its absence would lead to zero polar- 
ization: that is, leading collinear configurations are assumed for partons a and b 
inside unpolarized hadrons A and B, while transverse motion is considered in the 
fragmentation process. The final hadron, detected with momentum p A , is generated 
by the hadronization of a parton c whose momentum, p c = (p A — k±)/z, varies with 
k±; also the corresponding elementary process, ab — > cd, depends on k±. 

Pa is a function of the hyperon momentum p A = p L + p T and is normally 
measured in the AB cm. frame as a function of xf = 1"PlI \fs and px- 

Notice that, in principle, there might be another contribution to the polarization 
of a final hadron produced at large px in the high energy collision of two unpolarized 



hadrons; in analogy to Sivers' effect Jy], |TJ] one might introduce a new spin and 
dependent distribution function: 

A N f a y A (x a ,k ±a ) = f a y A (x a ,k ±a ) - f al/A (x a ,k ±a ) (10) 

= faT/A(x a , k ±a ) - f a y A (x a , -k_La) , 
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i.e. the difference between the density numbers f a t /Ai x a,k± a ) and f a i/ A (x a , k± a ) of 
partons a, with longitudinal momentum fraction x a , transverse momentum k± a and 
transverse polarization j or j, inside an unpo larized hadron A. 

This idea was first applied to unpolarized lepto-production |22| and to single 



spin asymmetries in pp' scattering pq|; the corresponding function, related to 



ftf /a( x o.i k± a ), was denoted by /i^. In the present case of transversely polar- 
ized A production this function would enter the cross-section accompanied by the 
transversity fragmentation function Hi(z) (or AD h u a r). We shall not consider this 
contribution here; not only because of the problems concerning A N f a y A (x ai k± a ) 
discussed below, but also because the experimental evidence that the hyperon po- 
larization is somewhat independent of the nature of the hadronic target suggests 
that the mechanism responsible for the polarization is in the hadronization pro- 
cessQ. A clean test of this should come from a measurement of Pa in unpolarized 
DIS processes, ip A T X gO . 



The k± dependent functions considered in Refs. []TTj, [T9|, [1^, ^3] (A N f a / A t, 
A N f a y A , A N D h / a i and A N D h y a , or, respectively, f^j,, hj;, and D^ T ) have the 
common feature that the transverse momentum direction is correlated with the 
direction of the transverse spin of either a quark or a hadron, via a sin <fr dependence, 
as can be seen from Eq. (0) for example. The reason for considering these functions 
is that this "handedness" of the transverse momentum compared to the transverse 
spin is displayed by the single spin asymmetry data in, for instance, pp^ — > n X 
and pp — > A^ X. However, the problem of using such functions is that naively 
they appear to be absent due to time reversal invariance. This conclusion would be 
valid if the hadronic state appearing in the definition of such functions is treated 
as a plane wave state. One can then show that the functions are odd under the 
application of time reversal invariance, whereas hermiticity requires them to be 
even. If, however, initial or final state interactions are present, then time reversal 
symmetry will not prevent the appearance of nonzero (naively) T-odd functions. In 
the case of a state like \Ph, Sh', X) final state interactions are certainly present and 
nonzero fragmentation functions A N D h / a t and A N D h i i a are expected. However, 
for distribution functions this issue poses severe problems and since we will only 



The main difference between the function A N D h / a r as originally proposed by 
Collins, and the function under present investigation A N D h r/ a , is that the former 



consider fragmentation functions here, we refer to Refs. [14], |22| for more detailed 
discussions on this topic. 

_ „x r. A/Vn 

'h/al 

is a so-called chiral-odd function, which means that it couples quarks with left- and 
right-handed chiralities, whereas the latter function is chiral-even. Since the pQCD 
interactions conserve chirality, chiral-odd functions must always be accompanied by 
a mass term or appear in pairs. Both options restrict the accessibility of such func- 
tions and make them harder to be determined separately. On the other hand, the 



1 This is not in contradiction with the observed spin transfer (Dnn) in pp^ — *• X ]2{|, which 
in the factorized approach can be described in terms of the ordinary transversity distribution and 



fragmentation functions, rather than in terms of A N D h 
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chiral-even fragmentation function can simply occur accompanied by the unpolar- 
ized (chiral-even) distribution functions, which are the best determined quantities, 
allowing for a much cleaner extraction of the fragmentation function itself. More- 
over, since chiral-even functions can appear in charged current mediated processes 



(as opposed to chiral-odd functions), more methods of extraction are available [[31 



As it was studied in Ref. the Collins function (or A N D h / a i) satisfies 
a sum rule arising from momentum conservation in the transverse directions. The 
same holds for the other fe^-odd, T-odd fragmentation function D^ T [53], 



EjdzJdkl^-Di T (z,k ± )=0, (11) 

or, in terms of the function A N D h y a , 

Y,C h M h = Y, fdzf d 2 k ± k ± sin<pA N D h]/a (z,k L ) = 0, (12) 
h h J J 

which is equivalent to Eq. flO) via Eq. (H). 

Notice that the above sum rule can be written as 

£ /W d 2 k±_ k ± D h y a (z, k±) = , (13) 

h 

and the same holds independently for D h y a (z, k±). Eq. ( |T3"1) has a clear nontriv- 
ial physical meaning: for each polarization direction (f or |) the total transverse 
momentum carried by spin 1/2 hadronsQ is zero. 

The sum over hadrons prevents a straightforward application of the sum rule to 
the case of A production alone. It can not be used as a constraint on the parameter- 
ization of the function to be fitted to the data. However, we note that the integral 
Ch Mh for each hadron type h is the same function of the energy scale (implicit in 
all expressions), apart from as yet unknown normalization. In this sense it closely 
resembles the tensor charge. In other words, the running of the functions are the 
same for any type of hadron and there is no mixing with other functions. The ratios 
for different types of hadrons are constants, which allow for checks of consistency 
between sets of data obtained at different energies, without the need for evolution. 
These constants are universal, if indeed the T-odd fragmentation functions are uni- 
versal. This universality is the main point of interest here: one wants to see if 
A polarization from unpolarized quark fragmentation is independent of the initial 
state, as it is implicit when writing down the factorized cross-sections Eqs. (§) and 
(H). At the present time, this can not be verified due to lack of data, but some 
predictions can be given [^] that will allow tests of this universality. 

We only consider the quark fragmenting into a A and not into other hyperons, 
like the S°. The latter actually produces a significant amount (20-30%) of the A's 

2 Strictly speaking, the sum over h is over all hadrons that carry transverse polarization, which 
might be true also for higher spin hadrons. 
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via the decay E° — > A7. The reason we do not introduce a separate S° fragmentation 
function at this stage is that the factorized approach by itself does not address such 
a separation (it is about a generic spin-1/2 hadron of type h), unless one introduces 
some additional input, like a model based on S77(3), or unless one applies it to 
separate sets of data for each hyperon (which are not available yet). Apart from 
that, for each quark flavour such a S° fragmentation function would evolve in the 
same way as the A fragmentation function, which implies that their relative fraction 
stays constant. In this way we can view the A fragmentation function as an effective 
fragmentation function that includes the contamination due to S°'s. Indeed, the 
fragmentation functions we shall use in next Section have been obtained from fits 
to inclusive A productions in e + e~ processes, independently of their origin. 

At a later stage one might make the distinction that the E° fragmentation is 
a different fraction of the effective A fragmentation function for different quark 
flavours. One can insert all this information with hindsight and correct for it, but the 
present approach cannot be used to acquire this information unless the data would 
clearly distinguish the A's coming from S°'s. Our approach of using an effective 
A fragmentation function would be more problematic if the S° would decay into 
other final states than A7 (which branching ratio happens to be 100%): then only 
a part of the total S° fragmentation function would be included into the effective A 
fragmentation function and this would be energy dependent. 

In the case of longitudinally polarized A production the S° — > A7 background 
gives rise to a depolarizing effect of about 10% pH , but in the present situation of 
transversely polarized A production this is not the case, since the photon does not 
carry away any of the transverse polarization and it hardly affects the definition of 
the plane compared to which the transverse polarization is measured. Therefore, 
the S° — > A7 background does not produce a significant depolarizing effect for the 
transverse A polarization and an effective A polarizing fragmentation function can 
be used also. 

We shall now consider both A and A production and attempt a determination 
of the polarizing fragmentation functions A N D A -\- i q by comparing results for Pa and 



P A from Eqs. (H) and (g) with data pS-p] 



2. Numerical fits of data on P A from pN — > AX processes 

Eq. (|8|), for proton-nucleon processes, can be rewritten as: 



E A d 3 aP N ^ AX 



rf3 PA a,b,c,d 

da ab ^ cd 



/ d 2 k± 




f dz j 




^min ** x 



x 2 b s 



X fa/p(Xa) fb/N^b) ~^ (x a ,X b ,k±) - {k ± -> -k±} 



A N D AU Az,k l 



(14) 



which deserves several comments and some explanation of notations. 

• In deriving Eq. (|T4| ) from Eq. (^) we have used the fact that A N D A i i c (z, k 
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Eq. (|3|), is an odd function of k±; the (+k±) integration region of k± runs over 
one half-plane of its components. 

The Xb integration has been performed by exploiting the 8(s + 1 + u) function; 
the resulting value of Xb is given by 

X a t$ t 

X b = ; — , (15) 

X a ZS + U® U 

where $ t and $ n are defined below. 

k± could have any direction in the plane perpendicular to p c ; however, due 
to parity conservation in the hadronization process - Eq. (|5[) - the only k± 
component contributing to the polarizing fragmentation function is that per- 
pendicular to P\, i.e. the component lying in the production plane, the x — z 
plane in our conventions. To simplify the kinematics we shall then consider 
only the leading contributions of k± vectors in the x — z plane. 

s, t and u are the Mandelstam variables for the pN — > AX process; in the 
simple planar configuration discussed above (p c and k^ both lying in the 
x — z production plane) they are related to the corresponding variables for the 
elementary processes by: 

s = 2p a -p b = x a x b s 

t = -2 Pa -p c =(x a /z)t$ t (±k ± ) (16) 
u = -2p b -pc= (x b /z)u<f> u (±k±) 



with 



* 9 t (±k x ) = g(k ± ) h T 2^±^ - [1 - 9(k±)\ ^ J (17) 
u * u (±k ± ) = g(k ± ) L ±2k ± ^ + [l- g(k ± )] ^ j (18) 



where g(k±) = J 1 — (k±/px) 2 and where ±/cj_ refers respectively to the con- 
figuration in which k^ points to the left or to the right of p c . At leading order 
in k±/pA one has: 



Pa V t pa V u 

The lower integration limits in Eq. fll4|) are given by: 



(19) 



u$ u (±k ± ) t$ t (±k±)+u$ u (±k ± ) 

ZS+t$ t (±k±) S 
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Notice that the integration limits are slightly different for +k± and — k±; when 
replacing k± with — k± inside the square bracket of Eq. (|14]) , one should not 
forget to change accordingly also the z and x a integration limits, and the value 
of Xb, Eq. (|15D, although the results are only marginally affected by this. 

Eq. ( |14D can be schematically written as 

da PN^X p A = daP N^X _ daP NM = £ U/pM ^ fb/N(xb) 

a,b,c,d 

<g> [da ab - cd (x ai x b , k ± ) - da ab ^ cd (x a , x b , -k ± )} ® A N D A r /c (z, k ± ) (21) 



which shows clearly that Pa is a higher twist effect, despite the fact that the 
polarizing fragmentation function A N D h i i a is a leading twist function: this is 
due to the difference in the square brackets, [da(+k±) — da(—k±)] ~ k^/p T , 
similarly to what happens for the single spin asymmetries in p^p — > irX El Q|. 



In the computation of the unpolarized cross-section E A d 3 a pN ~^ AX /d 3 p A in- 
trinsic transverse motion is significant only in limited phase space regions: we 
have checked that the values obtained in most of the kinematical regions of 
available data do not vary whether or not taking into account k±. Notice 
that when taking into account k± the expression for E\d 3 a pN ^ AX /d 3 p A is 



the same as Eq. (JA) with the — inside the square brackets replaced by a + 
and A N D A i/ c (z, k±) replaced by D\/ C (z, k±). 

When computing the cross-sections for scattering off nuclei, pA — > A^X, for 
which plenty of data are available, we have adopted the most simple incoherent 
sum, neglecting nuclear effects. That is, for the scattering off a nucleus with 
A nucleons and Z protons we use: 

da pA->AX = z da PP^AX + ( A _ Z ^ da pn^AX _ ^ 

We have checked that different ways of taking into account nuclear effects 
leave results for P\ - a. ratio of cross-sections - almost unchanged. The par- 
tonic distribution functions in a neutron are obtained from the usual proton 
distribution functions by applying isospin invariance. 



Eq. ( |l4l) holds for any spin 1/2 baryon; we shall use it also for A's, using charge 
conjugation invariance to obtain q — >• A fragmentation properties from q — ► A 
ones, which implies D^/q = D\/ q and A N D^ ^ = A N D^ i q . 



We now use Eq. QT3D in order to see whether or not it can reproduce the data 
and in order to obtain information on the new polarizing fragmentation functions. 
To do so we introduce a simple parameterization for these functions and fix the 



parameters by fitting the existing data on Pa and Pa [j3q] -|l39|. 
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We assume that A N D A y c (z, k±) is strongly peaked around an average value k ^ 
lying in the production plane, so that we can expect: 

/ d 2 k ± A N D A , /C (z, k ± ) F(k ± ) ~ A»D A y c (z, k° ± ) F(k° ± ) . (23) 

Consistently, since in this case F(k±) depends weakly on k^ when k± -C Pt, in the 
computation of the unpolarized cross-section we use: 

/ d 2 k ± D A/c (z, k ± ) F(k ± ) ~ \ D A/C (z) F(k° ± ) . (24) 

J (+k i ) 2 



The average k° ± value depends on z and we parameterize this dependence in a most 
natural way: 

Kz a (l-z) b , (25) 



M 

where M is a momentum scale (M = 1 GeV/c); in performing the fit we demand 

>, 2 = 



that K, a and b are constrained so that they satisfy the kinematical bound p 2 



(Pa ~ kD/z 2 — Pa-, fr° m which k\ < (1 — z 2 )p A and 

fc° (*) £ (PA)min v 7 !^ ^ (1 GeV/c) v 7 !^ , (26) 

which implies a > and 6 > 0.5. The values of K, a and b resulting from our best 
fit will have a clear physical meaning. 

We parameterize Aq D A r / c (z, kj_) in a similar simple form: we know that it has 
to be zero when k±_ = and z — 1; in addition, the positivity of the fragmentation 
functions - Eq. (|5p - requires, at any k± value, \A N D h u q (z, k±)\ < D h / q (z,k±). 
However, according to Eqs. (p3|), ( 24D and to take into account the mentioned [see 
comment after Eq. (|2"DD] difference of the integration regions (+k±) and (—k±), which 
is significant at the boundaries of the kinematical ranges (when px — k± and when 
p T ~ pTmax) we prefer to impose the more stringent bound \Aq D A r / c (z, k±)\ < 
D A / c (z)/2. Following Ref. ||4"0fl , this is automatically satisfied by taking: 



M 



z a "(l 



cjd d q «i(c q + d q y«+ d « 



D A/ q {z) 



N , K f^-^ RnM (27) 

9 4"d d q "/(c q + d q )^+ d » 2 



= N n Z^(l-zY q 



2 ' 

where we have used Eq. (J25|), and we require c q = a + a q > 0, d q = b + f3 q > 0, and 
|JV 9 '| K<1. 

We are now almost fully equipped to compute Pa and P A , let us briefly discuss 
the remaining quantities appearing in Eq. fll4]) . 
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We sum over all possible elementary interactions, computed at lowest order 
in pQCD. The polarizing fragmentation functions - parameterized as in Eq. 
([27j) - are supposed to be non vanishing only for A valence quarks, u, d and 
s; similarly for A valence antiquarks u, d and s. All contributions to the 
unpolarized fragmentation functions - from quarks, antiquarks and gluons - 
are taken into account. 



We adopt the unpolarized distribution functions of MRST |41[] . We have 
explicitly checked that a different choice makes no difference in our conclusions. 
We fix the QCD hard scale of distribution (and fragmentation) functions at 2 
(GeV/c) 2 , corresponding to an average value pt — 1.5 GeV/c. Since the range 
of pt values of the data is rather limited, no evolution effect would be visible 
anyway. 

We use the set of unpolarized fragmentation functions of Ref. |42j| , which 
allows a separate determination of D A / q and -Da/?) an d which includes A's 
both directly and indirectly produced; it also differentiates between the s quark 
contribution and the u and d ones: the non strange fragmentation functions 
D A / U = D A / d are suppressed by an SU(3) symmetry breaking factor A = 0.07 
as compared to D A / S . In our parameterization of Aq D A y q (z, k Q ± ), Eq. (|27|) , 
we use the same D A / q as given in Ref. [ 42H , keeping the same parameters a q 
and f3 q (c q and d q ) for all quark flavours, but allowing for different values of 
N u = N d and N s . 



We can now use pQCD dynamics, together with the chosen distribution and 
fragmentation functions, and the parameterized expressions of the polarizing frag- 
mentation functions, in Eq. (|14]) to compute Pa and P A , by comparing with data 
we obtain the best fit values of the parameters K, a, b, N u = N d , N s , c q and d q 
introduced in Eqs. (|25|) and (^). Notice that we remain with 7 free parameters, c q 
and d q being the same for all flavours. 

Our best fit results (x 2 /d.o.f. = 1.57), taking into account all data |35|-|38 
shown in Figs. 1-5. They correspond to the best fit parameter values: 



arc 



if = 0.69 a = 0.36 6 = 0.53 (28) 

N u = -4.30 N s = 1.13 c q = 6.58 d q = 0.67 . (29) 

Let us comment in some details on our results. 

In Fig. 1 and 2 we present our best fits to Pa as a function of px for different xp 
values, as indicated in the figures: the famous approximately flat p T dependence, 
for p T greater than 1 GeV/c, is well reproduced. Such a behaviour, as expected, 
does not continue indefinitely with px and we have explicitly checked that at larger 
values of px the values of Pa drop to zero: the shape of such a decrease, contrary 
to what happens in the pt region of the data shown here, strongly depends on the 
assumptions about the nuclear corrections. It may be interesting to note that this 



10 



fall-off has not yet been observed experimentally, but is expected to be first seen in 
the near-future BNL-RHIC data. 

Also the increase of \P\\ with xp at fixed pt values can be well described, as 
shown in Fig. 3; the two curves correspond to pt = 1.5 GeV/c (solid) and p? = 3 
GeV/c (dashed-dotted). 

The best fits of Figs. 1 and 2 are compared to p-Be data |3"5|-[]3"5|]; these are 



collected at two different energies, in the p-Be cm. frame, \/s ~ 82 GeV p6|-|38 



and y/s ~ 116 GeV Our calculations are performed at y/s = 80 GeV; we have 
explicitly checked that by varying the energy between 80 and 120 GeV, our results 
for Pa vary, in the kinematical range considered here, at most by 10%, in agreement 
with the observed energy independence of the data. 



Some data from p-p collisions are also available; in Ref. |$5| a linear fit to P\(xp) 
is performed, collecting all data with p T > 0.96 GeV/c, for an average value (p?) = 
1.1 GeV/c. In Fig. 4 we show these data and the linear fit (central line); the upper 
and lower lines show the fit error band. The solid line is our computation, at pr = 1.1 
GeV/c, with all parameters fixed as in Eqs. (|28"D and (f29|): our fit reproduces the 
data with good accuracy. 

In Fig. 5 we show our best fit results for Pa as a function of pr for different 



xf values, as indicated in the figure: in this case all data |36|, [38) are compatible 
with zero, with large errors, and the measured xf range is not as wide as for Pa, as 
expected from the lack of overlapping between A and nucleon valence quarks. 

The resulting values of the parameters, Eqs. fl28|) and (p9|), are very realistic; 
notice in particular that b essentially reaches its kinematical limit 0.5 and the whole 
function ( p5| ) giving the average k±_ value of a A inside a jet turns out to be very 
reasonable. 

Mostly u and d quarks contribute to Pa, resulting in a negative value of N u ; 
instead, u, d and s quarks all contribute significantly to Pa and opposite signs for 
N u and N s are found, inducing cancellations. 

In Fig. 6 we plot | A^ P A t / u ,d\ and A^ P a t i s as given by Eq. fl27|) with the best fit 
parameters (|28|) and (p9|) . We show, for comparison, also D A / U ^ d and D^/ s : notice 
how a tiny value of the polarizing fragmentation functions, in a limited large z 
region, is enough to allow a good description of the data. This also shows that 
taking into account only valence quark contributions to A^Pat / q - as we have done 
- is a justified assumption. 

A different set of unpolarized fragmentation functions can be found in the litera- 



ture |43]| : it holds for the quark fragmentation into A + A and gives P( A +A)/g rather 
than a separate expression of Da/, and Dx/ q ] it would be appropriate to compute 
the A + A single spin asymmetry: 



P 



da AT + da AT - da Ai - da Rl 



= (^a + SU)(m--7-t) • m) 
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However, since one knows from experiments on p — N reactions that in the kinemat- 
ical range of interest da A <C dcr A (and this is confirmed in our scheme, simply due 
to the dominance of q over q in a nucleon), one can safely assume 

^a+a Pa ■ (31) 

This set - differently from the one of Ref. - assumes SU (3) symmetry and takes 

Da/u = Da/6, = Da/s- 

We have also computed Pa with this second set of fragmentation functions; as 
in the previous case we have parameterized A N D A y q according to Eq. (0), with 
the same values of c q and d q for each flavour, but different values of N u>d and N s . 
We can equally well (x 2 /d.o.f. = 1.85) fit the data on Pa, obtaining fits almost 
indistinguishable from those of Figs. 1-4, with the best fit values: 

AT = 0.66 a = 0.37 6 = 0.50 (32) 

N u = -28.13 N s = 57.53 c g = 11.64 d q = 1.23 (33) 

Notice that, also in this case of SU(3) symmetric fragmentation functions D\/ q , 
and using only data on Pa, one reaches similar conclusions about the polarizing 
fragmentation functions A N D A r/ q '- N Uyd ^ N s and not only is there a difference in 
magnitude, but once more one finds negative values for A^ P a t i u d and positive ones 
for Aq D A y s . This seems to be a well established general trend. Plots analogous 
to those of Fig. 6, would show also in this case, Eqs. (pTD and (p3|), A^P a t/ s > 

l^o r -^AT/tt,dl- 

Very recently, new sets of quark and antiquark fragmentation functions into a A, 
based on a bag model calculation and a fit to e + e _ data, have been published [fW|]. 



Both a SU (3) flavour symmetric and a SU(3) symmetry breaking set are available, 
although at a rather too low energy scale (/i = 0.25 GeV). Nevertheless, we have 
also tried using these sets in our scheme to fit the data on Pa and P A . once more 
we can fit the data, better with the symmetric than the asymmetric set, and, again, 
negative A^ P a t i ud and positive A^ P a t i a are found, with A^ P a t u > \Aq P a t i ud \. 

3. Conclusions 

A phenomenological approach has been developed towards a consistent expla- 
nation and predictions of transverse single spin effects in processes with inclusively 
produced hadrons; we work in a kinematical region where pQCD and the factoriza- 
tion scheme can be used, but px is not much larger than intrinsic k±, so that higher 
twist contributions are still important. This applies to several processes for which 
data are available, like p^p — > irX []16 , and pN — > A^X. The single spin effect 



originates in the fragmentation process, either of a transversely polarized quark into 
an unpolarized hadron - Collins' effect |12|] - or of an unpolarized quark into a 
transversely polarized hadron - the polarizing fragmentation functions Single 



spin effects in quark distribution functions |TT| have also been discussed [14], [17 . 
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We have considered here the well known and longstanding problem of the po- 
larization of A hyperons, produced at large pr in the collision of two unpolarized 
hadrons. We have assumed a generalized factorization scheme - with the inclusion 
of intrinsic transverse motion - with pQCD dynamics; the new, spin and de- 
pendent, polarizing fragmentation functions A N D A t i q have been parameterized in a 
simple way and data on pBe — *> A^X, pBe — > A}X and pp — > A^X have been used 
to determine the values of the parameters which give a best fit to the experimental 
measurements. 

The data can be described with remarkable accuracy in all their features: the 
large negative values of the A polarization, the increase of its magnitude with xf, 
the puzzling flat pt > 1 GeV/c dependence and the a/s independence; data from p-p 
processes are in agreement with data from p-Be interactions and also the tiny or 
zero values of A polarization are well reproduced. The resulting functions A N D^ i q 
are very reasonable and realistic. 

Different sets of unpolarized fragmentation functions - either SU(3) symmetric 
or not - lead to very similar conclusions about these new polarizing fragmentation 
functions describing the hadronization process of an unpolarized quark into a po- 
larized A: they have opposite signs for u and d quarks as compared with s quarks 
and their magnitudes are larger for s quarks. They are sizeable - with respect to 
the unpolarized fragmentation functions - only in limited z regions, yet they can 
describe remarkably well the experiments. 

These polarizing fragmentation functions have a partonic interpretation and a 
formal definition, Eq. (0); they are free from the ambiguities related to initial state 
interactions which might affect analogous distribution functions and we expect them 
to be universal, process independent functions. Our parameterization of A N i q 
should allow us to give prediction for A polarization in other processes; a study of 
tp — > A 1 * X and e + e~ — ► A T X is in progress [ j30|| . 
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Fig. 1: Our best fit to Pa data from p-Be reactions, as a function of pt 
and for different xf bins, as indicated in the figure. Only some of the 
bins are shown; see Fig. 2 for complementary bins. The experimental 
results, |37]]-[|3!|, are collected at two different cm. energies, yfs ~ 82 
GeV and y^s ~ 116 GeV. For each a^-bin, the corresponding theoretical 
curve is evaluated at the mean xf value in the bin, and at ^fs = 80 GeV; 
the results change very little with the energy. See the text for further 
details. 
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Fig. 2: Our best fit to Pa data from p-Be reactions, as a function of pt 
and for different xf bins, as indicated in the figure. Only some of the 
bins are shown; see Fig. 1 for complementary bins. The experimental 
results, p7||-|]39||, are collected at two different cm. energies, \fs ~ 82 
GeV and \fs — 116 GeV. For each a^-bin, the corresponding theoretical 
curve is evaluated at the mean xf value in the bin, and at yfs = 80 GeV; 
the results change very little with the energy. See the text for further 
details. 
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Fig. 3: Pa data for p-Be reactions, as a function of xf and for different 
Pt bins, as indicated in the figure. The data are collected at two different 
cm. energies, y/s ~ 82 GeV and y/a ~ 116 GeV, H-|||. The two 
theoretical curves, evaluated at yfs = 80 GeV, correspond to pt = 1.5 
GeV/c (solid) and px = 3 GeV/c (dot-dashed). 
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Fig. 4: Experimental results for Pa in p-p reactions, as a function of 
xf, from Ref. [|36]]. All data with pt > 0.96 GeV/c are collected, and 
(Pt) = 1-1 GeV/c. Also shown is a linear fit to the data, taken from 
Ref. (central line); the upper and lower dot-dashed lines show the 
corresponding fit error band. The solid curve shows the theoretical com- 
putation, at pt = 1.1 GeV/c, with all parameters fixed as in Eqs. (28) 
and (p9|). 
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Fig. 5: Our best fit to data from p-Be reactions, as a function of pt 
and for different xp bins, as indicated in the figure. The experimental 
results |36|, |38| are collected at the cm. energies ~ 82 GeV. For each 
x i?-bin, the corresponding theoretical curve is evaluated at the mean x f 
value in the bin, and at ^/s = 80 GeV; the results change very little with 
the energy. 
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Fig. 6: Plot of |A^ AT/ J (= \A$D_ 



^AT/dl) and A$D Ar/s , as given by 
Eq. ( p7| ) with the best fit parameters ( p8| ) and (29). For comparison we 
also show the unpolarized fragmentation functions D^/ u (=-Da/<2) an d 

Da/s H3- 
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